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SOUND SCATTERING BY RIGID OBLATE SPHEROIDS, 

WITH IMPLICATION TO PRESSURE 
GRADIENT MICROPHONES 

Algirdas Maciulaitis,* John M. Seiner, 
and Thomas D. Norum 
Langley Research Center 

SUMMARY 

The objective of this investigation was to determine the frequency limit below which 
soimd scattering by a microphone body is sufficiently small to permit accurate pressure 
gradient measurements. To determine this limit, soimd pressure was measured at var- 
ious points on the surface of a rigid oblate spheroid illuminated by spherical waves gen- 
erated by a point source at a large distance from the spheroid. This insured an essen- 
tially plane sound field. The measurements were made with small pressure microphones 
flush mounted from the inside of the spheroid model. Numerical solutions were obtained 
for a variety of spheroidal shapes, including that of the experimental model. Very good 
agreement was achieved between the experimental and theoretical results. 

It was found that scattering effects are insignificant if the ratio of the major cir- 
cumference of the spheroid to the wavelength of the incident sound is less than about 0.7, 
this number being dependent upon the shape of the spheroid. This finding can be utilized 
in the design of pressure gradient microphones. 

INTRODUCTION 

It has long been recognized that effective and economical suppression of undesirable 
noise must proceed from an adequate knowledge of the mechanism and location of noise 
generation. Since many noise generation problems of current interest involve very com- 
plex, spatially distributed noise sources, new analytical and experimental techniques are 
being developed to provide the needed acoustic data. Essentially, these techniques require 
synthesis of the wave vector components of an acoustic field and measurement of its spa- 
tial scale. One such technique, which has been proven useful in measuring the flux of 
acoustic energy from a jet, requires accurate measurement of the cross covariance of 
the acoustic pressure gradient (ref. 1). 

Langley-Industry Research Associate, now at Grumman Aerospace Corporation, 
Bethpage, N.Y. 


There are several factors in the design of pressure gradient microphones that limit 
the frequency range in which accurate resolution of these pressure gradients can be 
obtained. Of these, consideration of the effect of scattering and diffraction by the micro- 
phone body is of prime importance and provides the motivation for the present investiga- 
tion. Previous work on this problem is very limited, is almost exclusively theoretical, 
and deals with scattering by pressure microphones. Noteworthy contributions in this 
category are those of Jones (ref. 2) on the semi-infinite rod of circular cross section 
and of Matsui (ref. 3), who solved the same problem as Jones for more realistic bound- 
ary conditions on the face of the microphone. No previous work dealing explicitly with 
the scattering problem in pressure gradient microphone applications has come to the 
authors’ attention. 

Experimental determination of the upper frequency limit imposed upon a pressure 
gradient microphone by scattering involves several difficulties which are primarily asso- 
ciated with measurement of small surface pressure differences. In addition, the experi- 
mental results are also subject to imperfections in generation of a plane wave field, devi- 
ation from free-field conditions, and scattering from extraneous objects, such as the 
model support. It is desirable, therefore, to be able to support the experimental results, 
as far as possible, with exact theoretical solutions within the framework of linear 
acoustics. 

The simplest configuration of a pressvire gradient microphone is a short circular 
cylinder. Unfortunately, an analytical solution for the scattering about this simple shape 
does not exist at the present time. The closest geometrical shape for which a solution 
does exist is the oblate spheroid. The similarity between these two body shapes improves 
as the cylinder length becomes very short, which is the case of interest for pressure gra- 
dient microphones. The reqviired analj^ical computation for the oblate spheroidal func- 
tions can be carried out with very high precision by utilizing computer programs recently 
developed at the U.S. Naval Research Laboratory (refs. 4 and 5). 

In the present investigation, surface pressure distributions on a rigid oblate spher- 
oid having major and minor axes of 25.40 cm (10.00 in.) and 8.08 cm (3.18 in.), respec- 
tively, have been measured for various sound incidence angles over the frequency range 
0.200 to 3.475 kHz. In this report, the results of these measurements are presented and 
compared with those analytically derived. Also presented are anal 3 rtical scattering pres- 
sure distributions for a zero-thickness disk and for an oblate spheroid twice as thick as 
the experimental model to show the effect of changing thickness. For the oblate spheroid 
used in the experiment, additional anal 3 d;ical results demonstrate the dependence of scat- 
tering on the distance between the sound source and the scatterer. Finally, scattering 
calculations for oblate spheroidal shapes ranging from a disk to a sphere are related to 
the problem of pressure gradient measurement by showing how the presence of the micro- 
phone will affect the pressure gradient of the sound field. 
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SYMBOLS 


Measurements and calculations were made in U.S. Customary Units. This paper 
presents physical quantities in the International System of Units (SI) with the equivalent 
values in U.S. Customary Units given parenthetically. 

A area, m2 (ft2) 

a one-half of the interfocal distance (spheroid); radius (sphere and disk), m (ft) 

B mn coefficient in series expansion. Pa (psi) 

F frequency parameter, defined by equation (22) 

g percent of projected area 

hn spherical Bessel function 

K shape calibration factor, defined by equation (12) 

k wavenumber, 2 tt divided by wavelength, m“l (ft“^) 

m summation index 

Nmn normalization coefficient, defined by equation (5) 

n summation, or exponential, index 

Pn Legendre polynomial 

p acoustic pressure. Pa (psi) 

Pq amplitude factor of incident pressure. Pa (psi) 

^mn oblate spheroidal radial function 

r distance between source and field points, m (ft) 

®mn oblate spheroidal angular function 
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oblate spheroid axis system 


X,Y,Z 

x,y,z Cartesian coordinates, m (ft) 

a incidence angle, deg 

Neumann factor, 1 when m = 0 and 2 otherwise 
77 angular oblate spheroidal coordinate 

Q CO -latitude angle in spherical coordinates, rad 

I radial oblate spheroidal coordinate 

n surface pressure augmentation due to scattering, p/pj 

a scattering parameter, defined by equation (14) 

(f) azimuthal oblate spheroidal coordinate, rad 

Subscripts: 


b back side of body 

f front side of body 

g fraction of projected area 

i incident wave 

0 source location 

s scattered wave 

1 body surface 
Superscripts: 

(1) first kind 
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third kind 


Prime denotes derivative with respect to variable. 
Bar denotes average. 


EXPERIMENTAL APPARATUS 

The experiment was conducted in the anechoic chamber of the National Bureau of 
Standards. The inside dimensions of the chamber are 10 by 6.7 by 6.7 m (33 by 22 by 
22 ft). The chamber is anechoic down to 50 Hz. The test consisted of measuring the sur- 
face pressures on an oblate spheroid exposed to the harmonic sovind field emanating from 
a 10-cm (4-in.) diameter midrange loudspeaker. The acoustic field produced by the loud- 
speaker was surveyed by a microphone traverse at the location where the oblate spheroid 
was later placed. It was found that over the frequency range of interest (0.200 to 3.745 
kHz), the free-field soimd pressure level across the 25.40-cm (10.00-in.) diameter of the 
spheroid varied by less than 0.5 dB. For the two lowest test frequencies (0.200 and 0.463 
kHz), the relative variation in sound pressure level was within 0.1 dB. The distance 
between the loudspeaker and the oblate spheroid was 3.50 m (11.5 ft). A view of the loud- 
speaker and oblate spheroid model moimted inside the anechoic chamber is shown in 
figure 1. 


Oblate Spheroid Model 

The cross section of the experimental model is an ellipse having a major axis of 
25.40 cm (10.00 in.), a minor axis of 8.08 cm (3.18 in.), and an interfocal distance 2a of 
24.09 cm (9.48 in.). This particular shape had been selected to provide a reasonable sim- 
ulation of a typical pressure gradient microphone having a thickness-to-diameter ratio of 
1:4. It had been anticipated from the work of Leitner (ref. 6) and Wiener (ref. 7) that for 
a given body shape, scattering would be strongly dependent on the ratio of a characteristic 
body dimension to the wavelength of the incident sound. Therefore, the advantages of the 
large model size were twofold: greater accuracy in measuring surface pressure distri- 
butions and the capability to obtain meaningful scattering data at frequencies within the 
range of the loudspeaker and the microphones. This made it possible to use small inex- 
pensive electret-type microphones mounted within the oblate spheroid. 

The model was fabricated in two parts from an aluminum alloy. A drawing of the 
oblate spheroid is shown in figure 2. The two parts of the model are held together by 
screws from the back side, which are not shown. To keep the surface as smooth as pos- 
sible, the internal microphone wells were connected to the model surface by holes 0.24 cm 
(3/32 in.) in diameter and 0.64 cm (1/4 in.) in depth. 
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Although 10 microphone wells were provided, only 1 to 7 were actually used in the 
experiment, and these cover a 90° angle from the axis of symmetry to the edge of the 
model. The other three had been originally intended for aiming the spheroid at the loud- 
speaker. However, this idea was later abandoned in favor of alinement with a laser beam. 
Because of space requirements, the two microphone holes near the thinnest part of the 
spheroid were not in the same azimuthal plane as the other five. Consequently, the data 
from these two points will not be presented other than for normal incidence. 

The oblate spheroid was suspended by a 1.5-m (5.0-ft) length of 1.3-cm (1/2-in.) 
pipe which could be turned aroxmd its axis or adjusted vertically by means of a simple 
threaded device (fig. 3). This arrangement provided the means by which measurements 
were obtained on the surface opposite the speaker. 

As shown in figure 1, 0.08-cm (1/32-in.) diameter guywires were used to insure 
proper alinement between the loudspeaker and the oblate spheroid. To accomplish this, 
a 6-mW He-Ne laser was installed in the loudspeaker holder and adjustments were made 
in the model suspension by observing the beam of light reflected from a small, front- 
surface mirror of laser quality located on the model. This mirror was adjusted prior to 
the test to be parallel to the diametral plane of the spheroid. In this way, the alinement 
accuracy of the experimental system was held within a minute of arc. 

Instrumentation 

Acoustic surface pressures were measured with miniature electret microphones. 

A general description of these microphones can be found in reference 8. The micro- 
phones were installed 0.64 cm (1/4 in.) below the model surface. The exterior dimen- 
sions of the microphone enclosure (which also housed a field effect transistor preampli- 
fier) were 0.71 by 0.48 by 0.41 cm (0.28 by 0.19 by 0.16 in.). The frequency response of 
these microphones was constant within ±0.5 dB between 0.200 and 3.500 kHz. The open- 
circuit sensitivity, in volts per microbar, varied somewhat from microphone to micro- 
phone, ranging between -56 and -62 dB. As discussed in the following section, each 
microphone was carefully calibrated for each test frequency under conditions simulating 
the geometry of installation within the oblate spheroid. 

A photograph of the inside of the spheroidal model is shown in figure 4. The micro- 
phones were held in the microphone wells by U-shaped mica retainer plates. In figure 4, 
the center microphone has been removed to expose the microphone hole and well. Micro- 
phones were installed with a 0.08-cm (1/32-in.) thick silicone rubber gasket sealing the 
periphery of the front face of the microphone. Sovind pressure levels in the absence of 
the oblate spheroid were measured with a 1.3-cm (1/2-in.) diameter condenser micro- 
phone, whose accuracy had been calibrated with a piston-phone calibrator. 
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Equipment used during the experiments to drive the speaker and to measure acoustic 
pressures is schematically shown in figure 5. Microphone outputs were passed through 
a 10-Hz bandwidth filter and then sequentially read off from precision rms voltmeters. 

Microphone Calibration 

All electret-type microphones used in the experiment were calibrated for both rel- 
ative and absolute sensitivity in two stages. In the first stage, relative pressure response 
over the working frequency range was measured by means of a cavity-type calibrator with 
a standard 1.3-cm (1/2-in.) condenser microphone as a reference. Microphones were 
mounted inside a special adapter which had been designed for insertion into the calibrator. 
The interior of the adapter was an exact facsimile of the microphone wells in the oblate 
spheroidal experimental model, including the 0.24-cm (3/32-in.) diameter by 0.64-cm 
(1/4-in.) long hole from the front of the microphone. This procedure established the 
relative sensitivity of each electret microphone. 

Several problems arose in establishing an absolute sensitivity for the electret 
microphones. Cavity-type calibrators are sensitive to the type of microphone inserted 
for calibration. Since the reference and electret microphone with adapters impress a 
different radiation impedance in the cavity, measurement of the absolute sensitivity is 
subject to a small error (on the order of 1 dB in sound pressure level). Therefore, the 
absolute sensitivity of the electret microphones was established by a second calibration 
stage. A solid oblate spheroid having dimensions identical with those of the experimental 
model was machined to hold a standard 0.64-cm (1/4-in.) condenser microphone in the 
location corresponding to the central electret microphone of the model. One electret 
microphone was then calibrated against the 0.64-cm condenser microphone in an anechoic 
room at a distance of 3.50 m (11.5 ft) between the loudspeaker and the model. By using 
the absolute sensitivity of this microphone in conjunction with the relative sensitivity of 
each electret microphone, the absolute sensitivities of all microphones were determined. 
This procedure insured the absolute sensitivity of the electret microphones within ±0.1 dB. 

ANALYTICAL FORMULATION 


Coordinate System 

The coordinate system relevant to the problem at hand is the oblate spheroidal sys- 
tem defined in terms of Cartesian coordinates by 


X = a(l + cos (j) 

y = a(l + sin (p ^ 

z = a^q j 


( 1 ) 
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with the ranges 
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In this system, the Z-axis is the axis of symmetry; the surfaces | = Constant, 

77 = Constant, and cj) = Constant are, respectively, confocal oblate spheroids of inter- 
focal distance 2a, confocal semihyperboloids of revolution, and semiplanes originating in 
the Z-axis, The limiting cases represent an infinitely thin disk (| = 0 , a = Disk radius) 
and a sphere (| — °°, a — 0 , a| = Sphere radius). In figure 6, a sketch of the oblate 
spheroidal system is given, and figure 7 illustrates how the experimental model was ori- 
ented in this coordinate system. 

The incidence angle ot shown in figure 6 is measured from the positive Z-axis. 
Throughout this investigation the incidence vector is specified to lie in the X-Z plane. 

The exact solutions presented in the following sections for scattering by an oblate 
spheroid and by a sphere can be found in reference 9 in terms of the velocity potential 
rather than the pressure. The assumed sinusoidal time variation is suppressed through- 
out the following discussion. 


Point Source Solution 

The incident pressure due to a point monopole source located at the position 
{^o>Vo>^o) is given by 


„ikr 

p. = p 

^o la- 


cs) 


where r is the distance between the source and the field point, and k is the wave num- 
ber of the sound. The resulting pressure on the rigid oblate spheroid is 


P = 


2 iPr 


kali + 1 11^=0 


00 00 

cos |^m(0 - 0o)j ^ 


^mn(”ii^^’i^o) ®mn(“ii^’^o) ®mn(“ii^»^) 


n=m 


Nmn(-ika) (-ika,i?i) 


( 4 ) 


where Rmn ®mn are oblate spheroidal radial and angular functions, respectively, 

and 


N mn (~ii^a) ~ y [^®mn("ii^>n 



8 


( 5 ) 



Plane Wave Incidence 


The incident pressure due to a plane wave traveling in the direction of the negative 
Z-axis is 


Pi = PoP'*^ 

The resulting pressure augmentation (n = p/Pi) on the surface of the rigid ohlate spheroid 
is given by 


n 2e^^ ^ i^ So^(-ika,-l) SQ^(-ika,t}) 

ka(l + n=0 No^(-ika) Rj,3)’(-ika,i^l) 

The corresponding result for a sphere with radius a is 

igikz V" ("i)'^(2n + 1)P (cos 0) 

n = ) ; — 2 (8) 

(ka)2 n=o h^^^ (ka) 

where 6 is the azimuthal angle in spheroidal coordinates and Pj^ and hj^ are the 
Legendre polynomial and spherical Bessel function, respectively. The derivation of equa- 
tion (8) can be foimd in reference 10. Since the results equivalent to equations (4) and (7) 
are presented in reference 9 without derivation, these derivations are included for the 
sake of completeness in the appendix. 


Effect of Scattering on the Pressure Gradient 

To determine the feasibility of using a given sized body as a pressure gradient 
microphone, one must determine the relationship between the pressure gradient in an 
acoustic field and the pressure acting on the body when placed within this field. A rea- 
sonable parameter to choose is the difference in pressure acting on opposite sides of the 
body divided by the body thickness. The question then arises as to how well this param- 
eter represents the pressure gradient that exists at the position of the median plane of the 
body when the body is removed from the field. 

The difficvilties involved in making this comparison can be divided into four distin- 
guishable areas. The first arises from the fact that a gradient can never be exactly cal- 
culated from a finite difference. Since this discrepancy exists independently of the pres- 
ence of the body, a term called the finite-difference factor is introduced to give an exact 
dependence of the pressure gradient on the free-field pressure difference. 
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The next two areas of difficulty concern differences associated with distortion of 
the field caused by scattering by the body surface. Even at low frequency, the pressure 
difference across the body will differ from that which exists over the same distance in a 
free field. This leads to the introduction of a so-called shape calibration constant which 
depends only on the shape of the body and which relates these two pressure differences. 
The second phenomenon due to scattering occurs at higher frequency and results in a 
deviation from the linear relationship between the frequency and the ratio of pressure dif- 
ference to pressure gradient. 

The final difficulty involves the quantization of the pressure difference across the 
active parts of the body. Even at low frequency, where the magnitude of the pressure is 
essentially uniform over the entire surface, curvature of the body results in a difference 
in phase at adjacent points on the body. Hence, an "average" pressure is introduced for 
both sides of the body. These pressures are determined by averaging the force in the 
gradient direction over the active portion of the body, i.e., that portion which would con- 
tain microphone sensors, and dividing by the area over which this force acts. The differ- 
ence in these average pressures will be dependent on the amount of surface considered to 
be active. However, as will be seen later, for a given body shape, this average pressure 
difference divided by the distance over which it occurs is independent of the amount of 
active surface at all frequencies for which the body would be considered for use as a pres- 
sure gradient microphone. 

To detail the discussion of pressure differences, consider a plane wave traveling in 
the negative z-direction, as shown in figure 8. The pressure difference between the planes 
z = Az/2 and z = -Az/2 is easily calculated from equation (6) to be 


Apj = -2ipQ sin 



Since the pressure gradient at the median plane, z = 0, is given by 


(9) 


= ikPo 

the relationship between the pressure difference and the pressure gradient can be 
expressed as 



( 10 ) 


_ sin (k Az/2) 

k Az/2 \9 z/2,_q 


The bracketed term represents the difference between the true pressure gradient and the 
finite pressure difference, and was previously termed the finite -difference factor. 
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When a body of thickness Az is placed in the field, the resulting pressure differ- 
ence across the body Ap will differ from Ap^. If a shape calibration constant K is 
defined as the ratio of these two differences at low frequency, 

/ApA 

K = lim — i (12 

k-0 \Ap j 

equation (11) can be rewritten as 


Az 


sin (k Az/2) 

k Az/2 _ \9z /^_Q 


(13) 


Although this equation is strictly valid only in the limit of zero frequency, deviations from 
it will occur only after a certain frequency is attained. Below this frequency, an essen- 
tially linear relationship is seen to exist between the pressure difference across the body 
and the pressure gradient in the free field. 

Deviations from this linearity arise at higher frequency because of two effects. 
First, as the frequency increases, the finite -difference factor begins to vary from its 
value of unity at low frequency. This merely expresses the fact that in order for a linear 
relationship to exist, the body thickness must be small in comparison with the wavelength 
of the sound, even if the body has no effect on the sound field. The second effect is due to 
high-frequency scattering, which causes the magnitude of the pressure on the body to be 
different from the magnitude of the free-field pressure. The frequency limit above which 
these effects become significant is investigated by introducing a scattering parameter ct: 


k Az/2 ^ 
sin (k Az/2) Az 



(14) 


This factor is chosen since by equation (13) it will be unity at low frequency. 

The value of the shape calibration constant and the behavior of the scattering param- 
eter can be determined for the oblate spheroid. For a plane soimd wave traveling in the 
direction normal to the major axis of the spheroid (fig. 9), the average pressure difference 
across the body is calculated in the manner previously indicated. For any arbitrary per- 
centage g of projected area, the average pressure acting on the front side of the spher- 
oid is 
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(15) 


Pf = 


C pdA 

^ 

r dA 
•Ja„ 


where dA is an element of projected area and Ag is the total area of interest. Since 
an element of projected area on the spheroid I = is given by 


dA = 27ra2 1 + )tj d?j 


this average pressure becomes 


(16) 



(17) 


where ? 7 g is the value of 77 


that encloses the area Ag, and is easily shown to be 



After the average pressure on the corresponding area on the back side of the spheroid ^ 
is computed in the same way, the pressure difference across the spheroid is calculated 
from 


Ap = Pf - Pjj 


(19) 


The body thickness corresponding to this difference in average pressure is computed in 
the same manner and is given by 


Az = 2 



z dA 


dA 


- "si 


( 20 ) 


From equations (10) and (19), the scattering parameter defined by equation (14) can 
be written 
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( 21 ) 


K .. I 

fPf - ^ 

2 sin (k Az/2) 

\ Po J 


For the oblate spheroid, the body pressure and thickness can be computed from equa- 
tions (17) and (20), respectively. 

Note that equation (21) can be used for two purposes. First, the shape calibration 
constant K can be determined from the fact that the scattering parameter a must 
equal unity at low frequency. Second, the frequency limit can be found by noting the vari- 
ation of the scattering parameter with a frequency parameter. Above this limit the pres- 
sure difference across the spheroid is no longer simply related to the pressvire gradient. 
The frequency parameter is conveniently chosen to be equal to the ratio of the body cir- 
cumference (normal to the wave) to its wavelength, which for the oblate spheroid is 


I 

F = ka^l + 


( 22 ) 


Computer Programs 

The oblate spheroidal fvinctions needed in the calculation of the pressure distribution 
on the spheroid were obtained from modifications of the computer programs described in 
references 4 and 5. The computed results were determined to an accuracy of at least five 
significant figures. Computations were performed on a Control Data (CDC) 6600 computer 
system. 


DISCUSSION OF RESULTS 

As a preliminary check of the validity of the computing program, a comparison was 
made with the results of Leitner (ref. 6) . This comparison is for a plane wave incident 
on a disk, which is the limiting case of an oblate spheroid. It is apparent from figure 10, 
where the pressure augmentations calculated from equation (7) are plotted, that the agree- 
ment with Leitner’s computations is excellent, except in the region near the edge of the 
disk. The small differences near the edge of the disk are probably due to the greater 
accuracy in the present investigation, which was made possible by the extremely accurate 
subroutines for oblate spheroidal functions taken from references 4 and 5. Also shown in 
the figure are Wiener’s (ref. 7) experimental results. These differ from the calculated 
results for two reasons: the finite thickness of the disk used in the experiment; and 
experimental errors, partly caused by an external probe microphone used to measure 
surface pressures. 
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Comparison of Experimental and Analytical Results 

The experiment was performed at frequencies from 0,200 to 3.475 kHz, correspond- 
ing to a range of frequency parameter F from 0.455 to 7.911. The model was turned 
both to change the incidence angle and, since all pressure taps were located on one face 
of the model, to measure the surface pressure on the back side. The experimentally 
determined surface pressure distributions are shown in figures 11 to 14. Also included 
in these figures are the theoretically predicted values, as computed from equation (4). 

The results presented are for five values of the frequency parameter F (0.455, 1.055, 
2.110, 4.219, and 7.911) and four incidence angles a (180°, 150°, 120°, and 90°). The 
ordinate II in these plots is termed the surface pressure augmentation due to scattering 
and is computed as the ratio of the rms value of the pressure at a given point on the sur- 
face of the oblate spheroid to the value at that point in the absence of the spheroid. With 
the exception of the data on the back side of the model at a = 150°, the agreement between 
experimental and analjdical results is seen to be very good. The reasons for the dis- 
agreement at a - 150° have not been found; unfortunately, it was not feasible to repeat 
this portion of the experiment. 

As seen in figures 11 to 14, only a small amoimt of scattering is present at the low- 
est test frequency (F = 0.455). At higher frequencies the change in the pressure distri- 
bution around the spheroid ranges from a very drastic one at a: = 180° to almost a uni- 
form small increase in pressure for a = 90°. It should be kept in mind that the pressure 
distributions shown are for one particular plane through the spheroid - only that normal 
to the X-axis. 

It is interesting to note that for a = 180° (fig. 11), both theory and experiment 
indicate a pressure augmentation greater than a factor of 2 on the front of the spheroid at 
F = 2.110. As the frequency is increased further, the peak values of E approach 2, as 
one would expect. Figures 11(d) and 11(e) show that the pressure can peak at locations off 
the center line and that more than one "loud spot" can occur on the back side. The latter 
finding, which is analogous to the diffraction ring pattern in optics, is not unexpected. 

Pressure Distributions in the X- Z Plane for Oblate Spheroids 

The experimental data were obtained essentially in a single plane normal to the 
X-axis (see fig. 7), which for arbitrary plane wave incidence provides only a partial pic- 
ture of the scattering effects. A more complete picture of these effects is given in fig- 
ure 15, which shows the calculated surface pressure augmentation in the X-Z plane at the 
same source positions and frequencies as in the experiments. 

These results were computed from equation (4) and are plotted in polar form against 
the argument cos"^ rj. For incidence at 180° and 90°, where the pressure distributions 
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are symmetrical, only one -half of the distribution is shown. It should be noted that the 
results for F = 0.455 (fig. 15(a)) are plotted on an enlarged scale. 

One can note from these figures that whenever a loud spot exists on the back side of 
the spheroid, it occurs almost exactly opposite the incident wave direction. However, the 
pressure peaks on the front side can be fovmd far removed from the incident direction. 
Figure 15(a) indicates that scattering does not appreciably alter the surface pressure 
distribution at low frequency. 

Variation of n With Distance Between Source and Model 

The computed results plotted in figures 16(a) and 16(b) show how the surface pres- 
sure augmentation would vary for the oblate spheroid model = 0.336j used in the exper- 
iment for different distances between the loudspeaker and the model. In these calculations 
the soxmd source was assumed to lie on the Z-axis. The pressure augmentations for a 
separation distance of 35.0 m (115 ft) are essentially identical with those for plane wave 
incidence. It is apparent from these graphs that the distance used in the experiment, 

3.50 m (11.5 ft), provided results very close to those for plane wave incidence. As the 
frequency increases, the small differences that do exist become progressively smaller, 
as can be observed by comparing figures 16(a) and 16(b). 

Effects of Doubling the Oblate Spheroid Thickness 

Calculated surface pressure augmentations for the experimental model and for a 
model having twice the thickness = 0.826j are compared in figures 17(a) and 17(b). 
These calculations are based on a sound source at o; = 180° and r = 35.0 m (115 ft). 

As expected, figure 17(a) demonstrates that the thickness of the body has little effect on 
the surface pressure distribution at low frequency; but at higher frequency (fig. 17(b)), 
the scattering effects become strongly dependent upon the geometrical shape. 

Effects of Scattering on Pressure Gradient Measurements 

One important consideration must be stressed before making inferences about the 
effects of scattering on pressure gradient microphones. All results obtained in the course 
of the present investigation, both experimentally and analytically, are for "hard" bodies, 
i.e., bodies with infinite acoustic impedance on their surface. Therefore, the findings 
presented are applicable either to microphones with very high impedance, such as the 
piezoelectric and piezoresistive t3q)es, or to those in which the active portions constitute 
a small fraction of the otherwise hard casing surface. 

By the use of equation (21) for low frequency where ct = 1, the shape calibration 
factor K was determined for the entire range of oblate spheroids from a disk to a 
sphere. These results are presented in figure 18 as a fimction of the ratio of minor to 
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major diameters of the spheroid. Note that K is always less than unity, which means 
that the pressure difference across the spheroid is always larger than the pressure dif- 
ference that would exist over the same distance in the absence of the spheroid. 

After the value of the shape calibration constant was determined, equation (21) was 
used to compute the variation of the scattering parameter a with frequency. For the 
particular oblate spheroid used in the experiment, the scattering parameter is plotted 
against the frequency parameter F in figure 19. Resxilts are shown for three choices 
of the percent of projected area over which the average pressure is computed. As can 
be observed, the scattering parameter is essentially independent of the percent of pro- 
jected area up to F = 1.0. A similar result was found for the sphere and other oblate 
spheroids. 

A useful feature of these results can be seen in figure 20. Here the scattering 
parameter is plotted against F for several oblate spheroidal shapes. As can be observed, 
the useful frequency range over which a body can be used to measure the free-field pres- 
sure gradient accurately is highly dependent upon the geometrical shape of the body. 

For example, in measurement of the pressure gradient within a relative error of ±1 dB 
(±12.2 percent), the results indicate that the high-frequency limit for a disk is approxi- 
mately only four -tenths of that for the oblate spheroid with | = 0.826. As shown in fig- 
ure 20, the disk represents the geometrical shape most sensitive to scattering effects. 

For this shape, incident wavelengths should be confined to values of F < 0.7 to obtain a 
relative error less than ±1 dB. 


CONCLUDING REMARKS 

Acoustic pressure distributions have been measured on the surface of a rigid 25.40- 
by 8.08-cm (10.00- by 3.18-in.) oblate spheroid in a harmonic plane wave field under 
anechoic conditions. Numerical solutions have been generated for the first time both to 
provide comparisons with experimental results and to investigate different oblate sphe- 
roidal shapes. Good agreement between experimental and analytical results has been 
obtained. 

The findings in this investigation indicate that pressure gradient microphones having 
oblate spheroidal shapes - all the way from a disk to a sphere - will be relatively free 
from scattering effects as long as the ratio of the major circumference of the body to the 
wavelength of the sound is kept below about 0.7. This information allows the selection of 
the proper size and shape for a spheroidal-shaped microphone in any intended frequency 
range. It was also found that below the frequency limit set by scattering considerations, 
the determination of the pressure gradient is dependent on the particular microphone 
geometry. This dependence must be either analytically or experimentally determined. 
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Since pressure gradient microphones are usually constructed in the shape of short 
circular cylinders, it is recommended that future scattering experiments using large 
models be extended to include such shapes. It is further suggested that in future work, 
small microphones be mounted flush with the model surface to avoid some of the calibra- 
tion difficulties encoxmtered in the present work. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
March 30, 1976 
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APPENDIX 


SOLUTION OF THE WAVE EQUATION IN TERMS OF 
OBLATE SPHEROIDAL FUNCTIONS 


Solution for Point Source 

The oblate spheroidal coordinate system allows an exact solution of the linearized 
wave equation. Although the solution for the pressure applicable to the scattering prob- 
lem for oblate spheroids follows a standard procedure and the final result is given in ref- 
erence 9, a step-by-step derivation could not be found in the literature and is therefore 
included herein for the sake of completeness. 

In brief, the solution involves the series expansion of the pressure of the incoming 
wave in terms of radial and angular oblate spheroidal fvinctions and trigonometric func- 
tions. Another series expansion containing undetermined coefficients is then assumed 
for the scattered wave. In this expansion the radial functions of the third kind are speci- 
fied, since they represent outgoing waves. The undetermined coefficients are evaluated 
by means of the boundary condition at the body surface 8p/9^ = 0 at ? = which 
states that the normal velocity is zero at the surface of the rigid spheroid. 

Suppression of the time variation allows the incident pressure due to a point source 
located at determined by equation (3). When the incoming wave is 

expanded in terms of the functions of interest, equation (3) becomes (see ref. 10) 


Pi ~ Z N — Hka ) ®nm("il^a,T?) Smn(-ika,T7o) ^^mn^"ika,i^) 
m n mn.' ' 


m,n 


X cos - 0 q) 


(Al) 


where is the radial oblate spheroidal function of the first kind and the other terms 

have been previously defined. 

Specifying the pressure due to the scattered wave as 




Bmn Smn(-ika,r/) 


m 


(0 - <Po) 


(A2) 


m,n 
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and using the boundary condition gives 


Bmn - 


-2l€mPo Ri^i,'(-lka,lit) 


3mn(“ika>Tjoj 


(A3) 


The pressure existing at any point on or outside the spheroid 1 = ^1 is then the 
sum of the incident and scattered pressures. On the surface of the spheroid, this addition 
yields 


2ip^ 


P = 


ka^l + 




m,n 


Smn(-ika,77) Smn{-^^,Vo) ^mn('^^>^^o) 


Nmn(-ika) 


(3)' 

mn 


) 


- cos 

m(<^ - <t>o) 

) 


(A4) 


where use of the Wronskian for the radial functions (ref. 10), which can be written as 


R 


( 1 ) 

mn 






1 


ka 1 + I 



has been used. 


(AS) 


Solution for Plane Wave Traveling in Negative z- Direction 

The incident pressure due to a plane wave traveling in the negative z-direction is 
given by equation (6). Expanding this in terms of the spheroidal functions yields (ref. 10) 


Y i SQj^(-ika,-l) SQj^(-ika,T 7 ) R^(-ika,i^) 


n=0 


(A6) 


The pressure due to the outgoing scattered wave is again expanded in terms of the radial 
functions of the third kind as 


Ps = Z ®0n 


n=0 


(A7) 
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The coefficients are evaluated as for the point source by the specification of zero pres- 
sure gradient normal to the body surface. Proceeding in the same manner as for the 
point source leads to the result 

p 2p^ ^ 5o„(-ika,-l) So^(-ika,77) 

ka(l + ^ 1 ^) n=0 No^(-ika) 

at all points on the surface of the spheroid ^ 
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Figure 6.- Oblate spheroidal coordinate system. 
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Figure 7.- Side view of experimental setup in terms of oblate spheroidal 
coordinate system for 180° incidence. 
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Figure 11.- Comparison of measured and calculated surface pressure augmentation due to scattering 

for 180® incidence to the oblate spheroid = 0.336. 
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Figure 12.- Comparison of measured and calculated surface pressure augmentation 
due to scattering for 150® incidence to the oblate spheroid = 0.336. 
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Figure 13.- Comparison of measured and calculated surface pressure augmentation 
due to scattering for 120® incidence to the oblate spheroid ^ < = 0.336. 




Figure 14.- Comparison of measured and calculated surface pressure augmentation 
due to scattering for 90° incidence to the oblate spheroid = 0.336. 
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Figure 20.- Variation of scattering factor with frequency parameter for several geometric shapes 
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